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ABSTRACT 

We introduce a new survey to map the radio continuum halos of a sample of 35 edge-on spiral 
galaxies at 1.5 GHz and 6 GHz in all polarization products. The survey is exploiting the new 
wide bandwidth capabilities of the Karl G. Jansky Very Large Array (i.e. the Expanded Very 
Large Array, or EVLA) in a variety of array configurations (B, C, and D) in order to compile 
the most comprehensive data set yet obtained for the study of radio halo properties. This is the 
first survey of radio halos to include all polarization products. 

In this first paper, we outline the scientific motivation of the survey, the specific science 
goals, and the expected improvements in noise levels and spatial coverage from the survey. Our 
goals include investigating the physical conditions and origin of halos, characterizing cosmic ray 
transport and wind speed, measuring Faraday rotation and mapping the magnetic field, probing 
the in-disk and extraplanar far-infrared - radio continuum relation, and reconciling non-thermal 
radio emission with high-energy gamma-ray models. The sample size allows us to search for 
correlations between radio halos and other properties, including environment, star formation 
rate, and the presence of AGNs. In a companion paper (Paper II) we outline the data reduction 
steps and present the first results of the survey for the galaxy, NGC 4631. 

Subject headings: ISM: bubbles - (ISM:) cosmic rays - ISM: magnetic fields - galaxies: individual 
(NGC 4631) - galaxies: magnetic fields - radio continuum: galaxies 
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1. Introduction 

We introduce a new survey to detect radio con- 
tinuum haloa^ in spiral galaxies using the Karl 
G. Jansky Very Large Array (hereafter, the Ex- 
panded Very Large Array, or EVLA). "Contin- 
uum Halos in Nearby Galaxies - an EVLA Sur- 
vey" (CHANG-ES) is observing 35 nearby edge-on 
galaxies (Table [TJ Fig. [!} in the radio continuum 
at two frequencies (1.5 GHz and 6 GHz, i.e. in 
L-band and C-band, respectively) in all polariza- 
tion products, and over a range of spatial scales. 
Our overarching goals are to characterize the na- 
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1 By 'halo', we mean the the gas, dust, cosmic ray, and mag- 
netic field components above a galaxy's plane, rather than 
stellar or dark matter halos. We will refer to the region, 
0.2 Si 2 ^ 1 kpc, as the disk-halo interface and use 'halo' 
for emission on larger scales (z 1 kpc). 'High-latitude' 
or 'extraplanar' are also used to describe either of these 
components. 



ture and prevalence of radio halos, their magnetic 
fields and the cosmic rays that illuminate those 
fields. By initiating a survey, we can address the 
connection between radio halos and the underly- 
ing galaxy disk as well as galaxy environment and 
provide a coherent data set with legacy value for 
the scientific community. 

In this first paper, our goals are to outline the 
motivation for the survey (Sect. [5]), to describe its 
scientific goals (Sect. [3]), to introduce the galaxy 
sample and its characteristics (Sect. [4}, and to de- 
tail the EVLA observations, including our specific 
goals for sensitivity, resolution and s patial cover- 
(Sect.[S]). In a companion paper (jlrwin et al. 



hereafter, Paper II) we present the first re- 
sults of the survey from test observations of the 
galaxy, NGC 4631. 

2. Motivation 

Since the early pioneering observations of edge- 

on galaxies dEkers fc Sancisi Hl977tlAllen. Baldwin fc Sancisi 



1978; Hummel fe van der Hulst lll989tlHummel. Beck, fe Dettmar 



19911 ). it has been clear that radio continuum 
emission can extend well above galaxy disks, that 
is, to z heights^ greater than 500 pc and often 
to many kpc. Radio continuum emission rep- 
resents just one component of the multiphase 
gaseous halos that are typically present. For 
example, halos have been observed in Ha emis- 
sion and other optical lines (called the 'extra- 
planar diffuse ionized gas' or eDIG), in soft X- 
ray emission, in the 21 cm HI line, in infra-red 



such as CO and PAHs (e.g. Bregman & Houck 


1997t iLee et al. 2001: 


Irwin & Seaquist 1990: 


Martin & Kern 


2001: 


Heesenetal. 1 2009allb: 


Whalev et al. 


2009: Wane et al. 19951 


2001; 


Howk & Savage 1 19991: Rossa & Dettmar 


2003; 


Rand 19961). The halo radio continuum emis- 



sion is primarily from synchrotro n radiation (e.g. 
Irwin. English fe Sorathial Il999l) implying the 



presence of cosmic rays (CRs) and magnetic fields 
out of the plane. 

Although radio continuum studies of edge- 
on galaxies have been conducted in the past 

(IHummel fc van der Hulst ll989HHummel Bec k, fc Dettmar 
1991 : Dumke et al. 1119951 llrwin. English fc Sorathia 



2 The height, z, is measured perpendicular to the plane of a 
galaxy disk. 
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19991: llrwin. Saikia. fc English I l2000h . the limited 
number of galaxies for some of the surveys, fre- 
quency coverage, sensitivity limits, and/or ability 
to detect a variety of spatial scales have made it 
difficult to systematically and consistently charac- 
terize the properties of radio halos. For example, 
those studies have observed (in the same order 
as referenced above) 4 galaxies at 1.5 GHz, 181 
galaxies at 5 GHz, 6 galaxies at 10.7 GHz, 16 
galaxies at 1.4 and 6 GHz, and 16 galaxies at 
1.4 GHz. The spatial scales that were detected as 
well as the sensitivity limits (see below) varied sig- 
nificantly. CHANG-ES is attempting to improve 
upon previous surveys in a number of ways. 

To begin with, a galaxy sample of significant 
size (Sect. [4j is being observed over a variety of 
angular spatial scales (Sect. [5]). The latter is par- 
ticularly important since previous observations of 
sufficient spatial resolution and sensitivity have 
shown that structure exists in halos on virtually 
any scale that is observed. The 'smooth' halos 
that are observed on broad scales resolve into dis- 
crete features on a variety of smaller scales. 

Low frequencies (1.5 and 6 GHz) have been cho- 
sen for their sensitivity to non-thermal emission 
and for the prospects of carrying out Faraday rota- 
tion analyses in galaxy disks and halos. Ours will 
be the first systematic survey for Faraday rotation 
in galaxy halos. The frequency choice also ensures 
that complementary data (for example, single dish 
data) can be obtained, where needed. 

CHANG-ES can also now exploit the wide 
bandwidths that are provided by the EVLA's 
Wideband Interferometric Digital ARchitecture 
(WIDAR) correlator. These wide bandwidths 
permit unprecedented improvements in contin- 
u um sensitivity. For example, the 5 GHz survey 



of iHummel. Beck, fc Dettmar 



(|199lh achieved a 
1 while the 1.4 GHz 



best rms noise of 80 ^J y beam 
and 5 GHz surveys ofE rwin. English fc Sorathia 



(Il999h achieved best rms values of 43 /iJy beam 
and 46 /xJy beam -1 , respectively. Our sensitivity 
limits should improve upon these previous surveys, 
in some cases by more than an order of magnitude 
(see Sect. [5]) . 

In addition, the wide bands provide a wealth of 
spectral information which will be unique to this 
survey. For example, it is now possible to measure 
a spectral index within a specific observing band 
(and possibly more than one for strong emission 



if there is curvature). We will therefore obtain at 
least three spectral indices from observations at 
two frequencies (Sect. [6]). 

In summary, CHANG-ES improves on previous 
surveys by observing a large number of galaxies 
over a variety of well-defined angular spatial scales 
at two frequencies, and to low noise limits. By ex- 
ploiting the wide-band capabilities of the EVLA, 
a wealth of spectral index information is available. 
Moreover, CHANG-ES is the first survey of edge- 
on galaxies for which all polarization products will 
be obtained, providing important information on 
halo magnetic fields. 

3. Science Goals 

3.1. Physical Conditions and Origin of 
Gaseous Halos 

The disk-halo interface and broader halo are 
dynamic entities that connect starkly different en- 
vironments: the normal ISM of a galaxy disk, and 

the rarefied intergalactic mediu m (IGM). A vari- 

ety of observational evidence (e.g. Dahlem. Lis enfeld. fc Golla 



19951: iRossa fc Dettmar ll2003l:IStrong et al. 112 004: 



Rossa et al. ll2008tlTullman et al. 1120061 ) points to 
a connection between the supernova rate per unit 
area in the disk and the intensity of the halo ra- 
dio continuum, X-ray and eDIG emission. In a 
few cases, it is also possible to associate discrete 
disk-halo features with underlying star forming 
regions, the latter also presumably being a source 
of eDIG ionization. 

Most models have therefore focussed on the 
critical link between extra-planar gas and star 
formation-related activity (largely supernovae) in 
the underlying disk. Such models have i n clude d 
galactic c himneys ( Norman fc Ikeuchi 1989) , 



fountains (IShapiro fc Field I 119761). s uperbubble 



blowout (|Mac Low fc Ferrara l Il999h and CR 



driven Parker instabilities (e.g.|B reitschwcrdt et al. 



1991 : Hanasz. Otmianowska-Mazur. fc Lesch Il2002t) 
CRs, in particular, appear to be a n essen- 



tial ingredient in driving outflows (e.g. Ipavich 
19751 : iBoulares fc Cox Ill99d : lEverett et aH l2008: 



^901 

Everett. Schiller fc Zweibel 1 12010| ) and it is this 
component, along with magnetic fields, that are 
critical to understanding outflow physics. In the 
disk, for example, CRs, the magnetic field, and 
hot gas have roughly equal pressures, whereas in 
the halo far from the plane, the CR and mag- 



3 



netic field pressures appear to dominate; however, 
it is not clear whether such statements are uni- 
versally true. Although we now have informa- 
tion on the strength and configuration of halo 
magneti c fields in so me galax i es [fo r reviews, see 



Krause 



(|2009h and iKrause I (|201l[ )]. an impor- 
tant goal of CHANG-ES is to consistently measure 
non-thermal pressures in galaxies with a variety of 
SFRs. Discrete disk-halo features can also be ex- 
amined for possible associations with star forming 
sites in the disk. 

Yet, in spite of the abundant evidence that links 
star formation to extra-planar gas, there are still 
some puzzles that suggest that underlying star 
formation is not a sufficient explanation for the 
existence of halos. There is evidence, for example, 
that a galaxy's global magnetic fiel d configuratio n 
seems not to be influenced by SFR ( Krause 1120091) . 
Optical line ratios suggest that a source of heating- 
other than photoionization from i n-disk sources is 



requi r ed to ionize the eDI G (e.g. ITullman et al. 



2000; Rand et al. 2011 ); possibilities include 



shocks, and heating from ma gnetic reconnection 
(jZimmer. Lesch. fc Birk 1119971 ) . In the latter case, 
for examp le, specific magnetic field topologies are 



predi cted (jHanasz. Otmiano wska-Mazur . fc Lesch 



2002). In another model (jGissinger. Fromang fc Dormv 
2009) vertical inflow plays a role. In still an- 
other, hydraulic jumps at the arm-interarm 
boundary produce l ong-lived vertical features 



(jMartos fc Coxlll998l) . SFR also does not seem 



to affect ra dio continuum scale heights ((Krause 
l2009U201lh . 

In addition, evidence for a direct relationship 
between HI supershells and SFRs is weaker. Ex- 
pansion energies of some HI supershells, where 
measurable, appear to be too high to be accounted 
for assuming reasona ble numbers of supernovae 
and stellar winds ([Spekkens. Irwin, fc Saikia I 



2004; ISilich et al. 1 120061 and references therein) 
and some HI supershells appear to be present 
outside of the st ar forming disks of galaxie s 



(jlrwin fc Seaauist I Il99d IChaves fc Irwin I l200ll) . 



prompti ng suggestions of oth e r sour ces of energy 
(but see IChakraborti fc Rav I ( 2011 ) ) or that the 
shells have been formed from im pacting external 



clouds (e.g lRand fc Stone Ifl996l) 



The fact that halos lag in comparison to the un- 



deriving disk (ISwaters et al. Il997 


: Tiillman et al. 1 


2000t Schaap. Sancisi. & Swaters 


2000tlLee et al. 1 



200ll : IFraternali. Oosterloo. fc Sancisi 112004 iHeald et ah 



2007), and the difficulty in explaining these lags 
via dynamical /ballisti c models with the gas orig- 
inating in the di sk ()Fraternali fc Binnev I [2006; 
Heald et al. 1120071 ) have led to interpretations that 
include the acc retion of external material onto 
galaxy disks (lOoster oo. Fraternali. fc Sancisi 



20071 : IFraternali et al. 1 120071 ). We need to un- 
derstand the effects of magnetic field p ressure and 
tensi on on the dynamics of cycling gas ( Benjamin! 
2002) in order to account for energies and velocity 
gradients in such systems. 

It is clear that circulation and/or outflow can 
have a dramatic effect on the evolution of the 
galaxy itself, from modifications of the metal- 
licity gradient to understanding the energy and 
mass budgets of galaxie s, to altered SFRs. Mass 



flux estimates thu s far (Brcgman & Pildis l fl99i 



Wang et al.~l Il995t IFraternali et al. 1 120021 ) imply 



that disk-halo outflows are responsible for moving 
a large amount of gas and may even result in extra - 



planar star formation (e.g. ITullman et al. 1 12003). 
That a disk-halo flow is responsible for some of 
the Galactic inte rmediate-velocity clouds (IVCs ) 
is also likely (e.g. IWakker fc van Woerden 1 19971) . 

Galactic 'feedback' is an important component 
in ACDM numerical models of galaxy formation. 
Whether feedback is in the form of central star- 
bursts, global starbursts, or active galacti c nuclei 
(AGN s) is not yet known (for a review, see lBaugh"! 
2006), but sometimes type la supernovae (objects 
that are normally associated with Population II 
objects within a galaxy), are implicated as an im- 
portant driver ( Li et al. 2011b : Wang 201dl and 
references therein). Observations of SO galaxies 
and bulges, for example, indicate that observed 
X-ray ha los are weaker th an feedback would sug- 
gest (e.g. Li et al. 2011bl ). Moreover some galax- 



ies with both weak AGN and star formation can 



also ha ve soft X-ray luminous halos (e.g. iLi et al. 
2011al ) . But the relations of such outflows to the 



stellar mass and to the galactic environment and 
evolution history remains uncertain. While radio 
continuum emission appears to be mainly related 
to Population I objects as described above, it is im- 
portant to verify whether this is always true and 
to search for deviations that may reveal a link to 
older components. 
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3.2. Cosmic Ray Transport and Wind 
Speed 

Theoretical and observational effort on the 
transport mechanism for cosmic rays (diffusive, 
convectiv e or both) has been carried out for some 
time fe.g.lStrong||l978HPohl fc Schlickeise71 ll99d: 
Duric. Irwin, fe Bloemen lll998T l. Cosmic rav elec- 
trons (CREs) can lose energy via a variety of 
processes, i.e. synchrotron, Inverse Compton, adi- 
abatic, ionization, non-thermal Bremsstrahlung, 
although the latter two can be neglected in nor- 
mal galaxies (jCondon I Il992f ) . The relative im- 



portance of the remainde r depend s on environ 



ment, but Hee sen et al. I (|2009al) have shown 



that the shortest timescales are associated with 
synchrotron and adiabatic losses in the outflow 
galaxy, NGC 253. An observational consequence 
(and one that CHANG-ES can look for) is a mea- 
surable steepening of the spectral index where 
adiabatic losses transition to radiative losses (see 
Sect. [6] for a discussion of the separation of ther- 
mal from non-thermal emission). An analysis of 
the variation of magnetic field strength and ra- 
dio spectral index with z height, (B(z) and a(z), 
respectively) can also provide a wealth of infor- 
mation on CR transport. 

Outflow velocity is an important datum for 
input to numerical models of galaxy formation. 
Since winds are thought to be driven by CRs and 
couple to the thermal outflow, a measurement of 
CRE bulk velocity should apply to the thermal 
wind speed as well. The CRE bulk speed can be 
determined from the radio scalehe ight at several 



freque ncies and the CRE lifetime (jHeesen et al. 
2009al ). Applying this method to more galax- 



ies and in more detail is an important goal of 
CHANG-ES. The variation of these quantities 
with position across the disk can also constrain 
outflow models and provide information on CR 
transport within the disk as well. 

Outflows from galaxies may also be the mecha- 
nism whereby the intergalactic medium becam e 



magnetized ( Kronberg. Lesch. fe Hopp 19991 ) 



The bulk CRE outflow speed in NGC 253, for ex- 
ample, achieves 300 km s _1 in some locations, i.e. 
greater than the escape velocity, implying that this 



outflow shoulc 


enrich 


(Heesen et al. 


2009a) 



Although the conditions 
in nearby galaxies may differ substantially from 



the high redshift universe, we can address very 
specific issues, such as how outflows scale with 
SFR, SFR/area, the presence of compact nucleii, 
and gravitational potential. 

3.3. Faraday Rotation and the Origin of 
Galactic Magnetic Fields 

A fundamental question which CHANG-ES 
may help to address is the origin of the field itself. 
Seed fields alone should wind up too quickly to 
account for the polarization angles that are ob- 
served in galaxies. The a — f2 dynamo is thought 
to be t he best alternative for generating observed 
fields (|Beck et al. Ill996l) . Although the mean- 



field dynamo might operate separately in the disk 
and halo regimes, the regime with the more active 
dynamo may 'enslave' the oth er, resulting in a 



coup l ing between disk a nd halo (|Moss fc Sokoloff 



20081: iMoss et al. Il2010l) 



Yet in the few cases where halo fields are 
detected, they can show geometries that are 
not easily explained by dynamos alone. In- 
deed, the observational evidence for regular (i.e. 
dynamo-generated) fields in galax y halos so far 
is rath er weak - possibly NGC253 (jHeesen et al. 
2009b) and indirect evidence from polarized in 
tensity asymmetries seen at 20 cm in 7 galaxies 
(|Braun. Heald fc Beck I lioioh . 7 araday rotation 
studies of our own Galaxy also show no coher- 
ent vertical field structure at the Sun's position 
(|Mao et al. Il2010t ). 



Rather, the observations appear to support 
dynamo-generated fiel ds in the disk t hat couple to 



outflowing winds (see Krause 20091) . For exam- 



ple, in NGC 4631, the field appears to be perpen- 
dicular to the disk where star formation is larger, 
but at larger radii where the SFR decli nes, th e 



field becomes parallel to the disk again (jKrause 



2004). 'X-shaped' magnetic structures that have 



been obse rved in edge-on galaxies have been repro 



duced by Hanasz. Woltanski. fc Kowalik (2009) 



in global MHD simulations of CR-driven dynamos. 
In addition, magnetic helicity conservation re- 
quires that outflows be present (IShukurov et al 



20061: ISur. Shukurov fc Subramanian 1120071) . The 
need for winds in explaining magnetic field struc- 
tures is certainly consistent with the abundant ev- 
idence for outflows that has been outlined above. 
Understanding the structure of the magnetic 
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field requires observations of Faraday rotation, 
which measures the strength of the ordered mag- 
netic field along the line of sight, and synchrotron 
polarization, which measures the o rdered mag 
netic field in the plane of the sky. ISoida et al. 



(|2011l ) provide an example of the application of 
rotation measure information to the study of an 
edge-on galaxy with a radio halo. By utiliz- 
ing many frequency channels over two frequency 
bands and obtaining polarization data over a vari- 
ety of spatial scales, CHANG-ES observations can 
disentangle these effects and (with models) even 
start to probe the 3-dimensional structure of the 
magnetic field. Rotation measure (RM) synthe- 



sis tec hniqu es described i n lBrentiens fc de Bruvn . 
(|2005h and iHeald et all (j2009j), which make use 



of the wide frequency response and many channels 
of the WIDAR correlator, are being employed. 

3.4. The FIR-Radio Continuum Correla- 
tion 

Young hot massive stars become supernovae 
which are believed to be the main sites of CR ac- 
celeration in galaxy disks, and hot massive stars 
are also important sources of photons that can 
heat the dust in disks giving rise to FIR emis- 
sion. This common origin is thought to be at the 
heart of the FIR-radio continuum relation which 
has been observed among and within gal axies (e.g. 
Helou et al."1ll985l:lMurphv et al. Il2006j) . 

Since the mean free path of a dust-heating 
photon is smaller than the diffusion length of a 
CR electron, a radio image of a galaxy will be 
a smoothed version of a FIR image. Retrieval 
of scale length information can provide important 
information on the star formation timescale and 



radiation intensity (see iMurphy et al. II2008T) and 



on th e degree of field ordering (jTabatabaei et al 
2012j). Moreover, by separating thermal from 



non-thermal radio emission (Sect. [5]) and correlat- 
ing the result with IR emission in different wave- 
bands, we can investigate the warm dust-thermal 
radio correlation and the cold dust-nonther mal ra- 

200.4 



dio correlation 



eg. 



Tabatabaei et al. 



Brar. Irwin fc: Saikia 



2007h . 



In halos far from sites of in-disk SF, one would 
expect the difference between mean free paths to 
be greater still, and that the FIR-radio relation 
may even break down altogether at some scale 
height. Together with IR images, CHANG-ES 



should be able to investigate to what extent the 
FIR-radio correlation is or is not valid in galaxy 
halos. Environment may also play an important 
part. The low density relativistic plasma compo- 
nent of galaxy halos can p rovide a sensitiv e trace r 

~ rf2009h . 



Murphy et al. 



to environmental effects, 
for example, have shown that the relation is mod- 
ified in Virgo Cluster spirals, likely because of 
galaxy motions with respect to the intra-cluster 
medium. The CHANG-ES sample includes galax- 
ies within a wide range of environments. 

3.5. CRs and High Energy Modeling 

Synchrotron emission results from the electron 
component of CRs and so there is an tight con- 
nection between the gamma-ray emission that is 
produced by CR nuclei and lepton interactions 
with gas and radiation fields, and the radio con- 
tinuum emission that CH ANG-ES will observe . 
The G ALPROP codM dMoskalenko fc Strong 



1998; Moskalenko et al 



2002HPorter et al. Il2008t 



Strong et al. ll200dl20oil2009Ll2010h self-consistentlv 
links these and other components in a model of the 
Galaxy, predicting (polarized and total) broad- 
band fluxes from radio to gamma-ray energies; it 
also allows for the inference of other parameters 
such as the escape probability of particles and the 
diffusion coefficient. 

With modifications, this model will be ap- 
plied to other galaxies to predict synchrotron 
scale heights together with gamma-ray intensi- 
ties from each galaxy. Gamma rays give inde- 
pendent information on the energy density of the 
CRs and hence allow the determination of mag- 
netic field strengths without the assumption of 
energy equipartition. With the ability to put 
strict observational limits on galaxy radio halo 
sizes (one of the largest sources of uncertainty 
in GALPROP models of the Milky Way) we will 
be able to constrain the CR diffusion and energy 
losses for all species (primary /secondary electrons 
and nuclei). At present, NGC 253 and M82 are 
the only edge-on galaxies that have been detected 
in high-energy gamma rays (Acero et al. 20091: 



Abdo et al.H2010t lAcciari et al. I l2009h . but our 



gamma-ray flux predictions can be constrained by 
current upper limits. 

X-ray emission generated from inverse-Compton 



3 http:/ /galprop.stanford.edu/ 



G 



scattering also provides important relevant data. 
For example, existing X-ray observations from 
Chandra and XMM-Newton may already be use- 
ful for measuring or tightly constraining the hard 
X-ray (e.g., 2-10 keV) intensity and its spatial dis- 
tributions in galactic halos. The combination of 
the X-ray and the radio measurements can be used 
to measure the magnetic field and CRe density. 
Future X-ray satellite^ may provide additional 
constraints. 

3.6. Disks and Nuclei 

Radio continuum emission is a tracer of massive 
star formation and does not suffer from obscura- 
tion as do many other tracers such as Ha emission. 
This is particularly important for galaxies that are 
edge-on to the line of sight and so have high optical 
depths at other frequencies. For example, the Ha 
SFR, corrected for ob scuration using A24 /im data 
( Calzetti et al. has not been well-calibrated 

for edge-on galaxies. With correctio ns for contri- 



butio ns from lower mass stars (e.g. Ilrwin et al 



20111 ). the global radio emission in the disks of 
the CHANG-ES galaxies should therefore provide 
a good measure of total SFR which can be com- 
pared to SFRs from other tracers. The radio spec- 
tral index maps can also provide information on lo- 
cations of star formation and CR transport within 
the disk (see Sect. 13. 2[) . 

Compact radio cores are indicative of the pres- 
ence of AGNs and we also know that some spi- 
ral ga laxies exhibit nuclear jets and/or o u tflows 
(e.g. Hummel, van Gorkom. fc Kotanvil J l983: 
~Ulvestad 1119831 Ilrwin fc Seaauist Ill988l ). 
ulges are believed to be correlated 



IWilson 



with central black holes ( 


Gcbhardt et al. 2000: 


iFerrarese & Merritt 2000l) but buleeless svstems 


can also harbour AGNs ( 


Satvapal et al. 2009). 



Although we have some indication as to whether 
nuclear activity is present in our galaxy sample 
(see Table Q] 'Nuclear Type'), the presence of 
aligned emission or compact radio cores should 
help to distinguish whether nuclear activity is 
present. These results can be compared to in- 
formation at other wavebands and can also be 
correlated with galaxy environment. 



4 The Chinese Hard X-ray Modulation Telescope, NASA's 
Nuclear Spectroscopic Telescope Array (NuSTAR), and the 
Indian X-ray satellite, ASTROSAT, are possible examples. 



4. The Galaxy Sample 

The sample was tak en from the N earby Galax- 
ies Catalog (NBGC, iTullv I Il988l) which lists 
galaxy parameters in a consistent fashion and also 
provides other useful data such as a parameter 
indicating local galaxy density. The criteria for 
including a galaxy in our sample were: 

a) Inclination, i > 75°, to ensure a nearly edge-on 
aspect, 

b) Declinations, 5 > —23°, for accessibility (with 
adequate uv coverage) from the EVLA, 

c) Blue isophotal diameters, 4 < d25 (arcmin) < 
15; the lower limit ensures sufficient spatial res- 
olution and the upper limit avoids large galaxies 
that would require many pointings and mosaicing. 

d) Flux densities at 1.4 GHz, Si. 4 > 23 mJy, to 
ensure the likelihood of a detection. 

In addition, 3 galaxies which fell just outside 
of the above limits were also added because they 
were known to have extraplanar emission and for 
which good ancillary data were available. These 
were NGC 5775 (d 25 < 4 arcmin), NGC 4565 
(d 25 > 15 arcmin), and NGC 4244 (d 25 > 15 
arcmin and Si. 4 < 23mJy). There are 35 galaxies 
in total. 

Table Q] illustrates the range of properties of the 
galaxies. The distances range from 3.3 to 42 Mpc, 
the galaxies' morphological types span the range 
of spirals, both barred and unbarred and a variety 
of nuclear activity types are represented. They are 
also in a variety of environments (see p in Table[l}, 
from apparently isolated galaxies (e.g. NGC 2683) 
to Virgo Cluster members. A very nice aspect of 
this survey is that, other than the 3 galaxies men- 
tioned above, we have not specifically targetted 
galaxies with known halos and, as a result, there 
are some galaxies in the sample which have not 
previously been searched for halo emission. 

The choice of systems with detectable radio 
flux suggests that our galaxies are actively form- 
ing stars, although contributions from radio cores 
are also present in some cases. Star formation 
rates (SFRs, Table [T]), however, are modest, the 
maximum being 20 M yr^ 1 (NGC 4666) or, 
among those with exclusively HII region nuclear 
spectra (i.e. no known compact nuclear activity) 
the maximum is 6.4 M yr" 1 (NGC 5792). By 
comparis on, a sample of 417 nearby sta r forming 
galaxies (Moustaka s fc Kennicutt ll2006l ) has SFRs 
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ranging from 0.03 to approximately 170 M Q yr 



(jKennicutt 1 12008). Thus, our choice of galaxies 



with radio continuum fluxes above 23 mJy has not 
imposed a bias to high SFR in our sample. We are 
probing the halos of 'normal' galaxies. 

5. EVLA Observations 

CHANG-ES has been awarded up to 405 hours 
of observing time in dynamic scheduling mode at 
the EVLA. Observations are being carried out 
at 1.5 GHz and 6 GHz with 500 MHz and 2 
GHz bandwidths, respectiveljU. EVLA contin- 
uum observations are now made in spectral line 
mode which facilitates radio frequency interfer- 
ence (RFI) mitigation and also allows for multi- 
frequency synthesis techniques to be applied dur- 
ing imaging (see Paper II). We are using 2048 spec- 
tral channels at 1.5 GHz and 1024 channels at 6 
GH jf| which also permits the effective application 
of RM synthesis techniques (Sect. 133)) . 

The distribution of time over frequencies and 
arrays is given in Table [2j These values are best- 
case estimates, given that the observing time has 
been awarded in 'Shared Risk' moduli during the 
extensive commissioning period of the EVLA. All 
polarization products (Stokes I, Q, U, and V) are 
being obtained. 

Observations in the B, C and D array configu- 
rations at 1.5 GHz and in C and D configurations 
at 6 GHz ensure that a variety of spatial scales 
are probed, as outlined in Sect. Table [5] shows 
that the linear resolution spans a range from 56 pc 
to 9.4 kpc, depending on galaxy distance, observ- 
ing frequency and array configuration. In fact, the 
new wide bandwidths provide even more flexibil- 
ity than the table implies since the resolution will 
vary between the upper and lower ends of each 
observing band. With appropriate uv weightings, 
the array can also be 'tuned' to a wider range of 
spatial resolutions. 

The highest linear resolutions range from 69 pc 
to 0.88 kpc at 1.5 GHz and 56 pc to 0.71 kpc at 



5 Radio frequency interference may reduce these values some- 
what in practice (see Paper II). 

6 Hanning smoothing will effectively set the spectral resolu- 
tion at twice the channel width. 

7 Time has been granted under the 'Resi- 
dent Shared Risk Observing' (RSRO) program 
(https : //science .nrao . edu/f acilities/evla/early-science/r 



6 GHz, depending on galaxy distance. For sci- 
ence which requires a comparison between galax- 
ies at similar linear resolutions, we will either spa- 
tially smooth the data, apply appropriate mod- 
els, and/or limit the number of comparison galax- 
ies, as necessary. An example is a comparison of 
vertical scale heights between galaxies. A spatial 
resolution of approximately 500 pc should permit 
the measurement of a scale height, assuming an 
exponential decline. At 1.5 GHz, this is possi- 
ble for all galaxies closer than 24 Mpc (25 out of 
35 galaxies) and at 6 GHz, galaxies closer than 
30 Mpc (30 out of 35 galaxies) meet this crite- 
rion. Note, however, that we may detect emission 
to much g reater distances; for exampl e, the A 20 
cm map of Hummel fc Dettmar ( 1990f) shows the 
spectacular radio continuum halo of NGC 4631 be- 
ing traced to at least 10 kpc from the plane. 

On the other hand, the largest angular size de- 
tectable by the array is fixed by the shortest an- 
tenna spacing. At 1.5 GHz, this is 16 arcmin which 
is comparable to the optical diameter of our largest 
galaxies (Tabled]), but at 6 GHz, the largest de- 
tectable angular size is only 4 arcmin. Therefore, 
missing flux must be supplied by supplementary 
single-dish observations of the largest galaxies at 
1.5 GHz and of all galaxies at 6 GHqf|. This is 
especially important in order to derive accurate 
spectral index maps between 1.5 GHz and 6 GHz 
and also to derive proper scaleheights at 6 GHz. 
Details as to how the single dish data will be com- 
bined with the EVLA data will be described in a 
future paper. 

At 1.5 GHz, the full width at half maximum 
of the primary beam (Ops) is 30 arcmin so only 
one pointing is necessary for all galaxies. At 6 
GHz, however, Opb — 7.5 arcmin. Considering 
the sensitivity of the array to emission outside of 
the primary beam as well as observing time re- 
quirements, we have opted to observe all galaxies 
with d25 > l.39ps with two pointings. These 
pointings are placed on the major axis on either 
side of the nucleus, separated by one-half of the 
primary beam. 

The noise limits (Table [2]) improve over previ- 
ous observations at comparable spatial resolution. 



8 At the time of writing, we anticipate that these data may 
come from the Greenbank Telescope, which will have wide 
^..bandwidth receivers installed in 2012. 
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Note that the estimated confusion limit applies to 
the total intensity image and will be lower for Q 
and U images. In addition, it should be possible to 
make corrections for confusion in the low resolu- 
tion data (e.g. D array) using the high resolution 
(B array) results. 

A description of the data reduction procedures 
for EVLA data are provided in Paper II. 

6. On the Separation of Thermal and Non- 
thermal Emission 

Globally, we know that radio continuum emis- 
sion in galaxies is dominated by the non-thermal 
component at our observing frequencies. However, 
the thermal component is not negligible; for exam- 
ple, thermal contributions are typically 12% and 
26% of the total flux at 1 .5 GHz and 6 GHz, re- 
spectively (jCondon II1992I ). Since the bulk of the 
global emission is from the galaxy's disk rather 
than the halo, these fractions should also be typ- 
ical of disks. Few measurements of spectral index 
have so far been made in fainter galaxy halos; how- 
ever, in NGC 4631, halo em ission shows a steeper 



spect ral index than the disk (jHummel &: Dettmar 



1990), implying that the thermal contribution is 
weaker and possibly negligible in halos. 

The highest linear spatial resolutions for CHANG- 
ES data (~ 60 pc, Table[2]) are sufficient to resolve 
larger HII regions in normal galaxies (for exam- 
ple, the largest HII regions in M 31 e xceed 100 pc 
Azimlu. Marciniak. fc Barmbv I l201ll) . Therefore, 
the thermal fraction could be significant and, in- 
deed, even dominant in some beams in the disk 
and possibly the disk-halo interface. In Paper II, 
for example, we show that the simplest explana- 
tion for the spectral curvature that is observed in 
the disk of NGC 4631 at 1.5 GHz is a contribu- 
tion from thermal emission; we also show that two 
star forming complexes have flat spectral indices, 
consistent with thermal emission. 

Wideband EVLA observations (together with 
single-dish data to fill in missing spacings) con- 
tain a wealth of information that can help with 
the thermal/non-thermal problem. With multi- 
channel observations in two frequency bands, it 
is possible to obtain at least 3 spectral index 
measurements, at least one within each band, 
and one between the two bands at every point. 
Since a thermal spectrum is characterized by a 



known spectral index (I„ cx v~ 0A ) which does 
not vary with frequency, a measurement of 3 spec- 
tral indices may be sufficient to separate thermal 
from non-thermal emission in the case of a con- 
stant non-thermal spectral index, a^T- New al- 
gorithms, however, are permitting the measure- 
ment of changing spectral indices within observing 
bands (see Paper II for details) . For positions with 
sufficient signal-to-noise, curvature in olnt can be 
detected. This is particularly important for mod- 
elling cosmic ray propagation (e.g. Sect. 13.21) . 

7. Summary 

In this paper, we have introduced a new sur- 
vey, called CHANG-ES, to observe radio contin- 
uum halos in 35 edge-on, normal spiral galaxies. 
The galaxies are being observed in all polarization 
products in two different bands, 1.5 GHz (L-band) 
and 6 GHz (C-band), and over 3 different EVLA 
array configurations. This is the first comprehen- 
sive radio continuum survey of halos to include all 
polarization products. 

Our scientific goals include understanding the 
physical conditions and origin of gaseous halos, 
probing cosmic ray transport, and determining 
outflow speeds and the ultimate fate of the gas. 
Mapping the halo magnetic field configuration and 
investigating how fields couple to outflows are par- 
ticularly important goals, as is examining the FIR- 
radio continuum relation in halos. We will also 
measure synchrotron scale heights and attempt to 
reconcile these with models that focus on high en- 
ergy (including gamma-ray) emission. Ancillary 
data, such as infra-red and X-ray data will be ap- 
plied to some of these questions. 

This survey exploits the new capabilities of the 
EVLA, especially the wide bandwidths that con- 
tribute to reducing on-source integration times in 
the detection of galaxy halos. With up to 405 
hours of observing time, our combination of num- 
ber of galaxies observed, range of spatial scales 
probed, and low noise levels (see Tables Q] and 
[2]) significantly improve upon previous surveys of 
edge-on galaxies carried out at these two frequen- 
cies; in the case of noise levels, the improvement is 
more than an order of magnitude. Spectral index 
variations can also be observed and used to sepa- 
rate thermal from non-thermal emission as well as 
to model non-thermal spectral index variations in 
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regions of sufficiently high signal-to-noise. 

Our first CHANG-ES test observations of 
NGC 4631 at C array, along with details as to 
the data reduction procedures, are reported in 
Paper II. 
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Fig. 1. — The CHANG-ES logo, showing optical images of the galaxy sample and the 'Y' shape of the EVLA. 
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Table 1 

The CHANG-ES Galaxy Sample 



Name 


RA (J2000) a 


DEC (J2000) a 


i b 


d 2 5 c 




D e 


Typo f 


Nuclear 


Si. 4 ( 


ref) g 




SFR' 


io g (M T y 


P k 




(h m s) 


(d m s) 


(deg) 


(') 


(km/s) 


(Mpc) 




Type a 


(mjy) 


(10 9 L Q ) 


(M yr" 1 ) 


(log(M )) 


(Mpc -3 ) 


N 660 


01h43m02.40s 


+ 13d38m42.2s 


77 


7.2 


856 


12.3 


SBa 


HILLINER 


373 


(A) 


16.8 


5.7 


10.73 


0.12 


N 891 


02h22m33.41s 


+42d20m56.9s 


84 


12.2 


529 


9.48 


Sb 


HII 


286 


(C) 


12.1 


4.7 


11.19 


0.55 


N 2613 


08h33m22.84s 


-22d58m25.2s 


85 


6.8 


1679 


23.4 


Sl> 


HII 


69 ( 


B) 


9.68 


4.6 


11.58 


0.15 


N 2683 


08h52m41.35s 


+33d25ml8.5s 


7!) 


9.1 


415 


6.27 


Sb 


LINER, Sy2 


84 ( 


B) 


0.86 


0.36 


10.81 


0.09 


N 2820 


09h21m45.58s 


+64dl5m28.6s 


90 


4.1 


1576 


26.5 


SBc 




46.6 


(A) 


5.91 


2.0 


10.87 


0.28 


N 2992 


09h45m42.00s 


-14dl9m35.0s 


90 


4.4 


2212 


34.0 


Sa 


Syl,Sy2 


226 


(A) 


16.6 


7.2 


11.24 


0.21 


N 3003 


09h48m36.05s 


+33d25ml7.4s 


90 


6 


1480 


25.4 


Sbc 




30.4 


(A) 


3.87 


1.5 


10.75 


0.64 


N 3044 


09h53m40.88s 


+01d34m46.7s 


90 


4.4 


1335 


20.3 


SBc 


HII 


103 


(A) 


7.14 


2.6 


10.81 


0.19 


N 3079 


10h01m57.80s 


+55d40m47.3s 


88 


7.7 


1125 


20.6 


SBcd 


LINER, Sy2 


821 


(C) 


39.1 


13 


11.27 


0.29 


N 3432 


10h52m31.13s 


+36d37m07.6s 


82 


4.9 


611 


9.42 


SBm 


LINER, HII 


120 


B) 


1.31 


0.39 


10.19 


0.19 


N 3448 


10h54m39.24s 


+54dl8ml8.8s 


78 


4.9 


1391 


24.5 


SO-a 




51.3 


(A) 


6.65 


2.1 


10.95 


0.24 


N 3556 


llhllm30.97s 


+55d40m26.8s 


81 


7.8 


697 


14.09 


SBc 


HII 


217 


(A) 


12.5 


4.9 


10.78 


0.15 


N 3628 


Ilh20ml7.01s 


+ 13d35m22.9s 


87 


14.8 


846 


8.50 


Sb 


HILLINER 


470 


(C) 


6.99 


2.3 


11.17 


0.39 


N 3735 


Ilh35m57.30s 


+70d32m08.1s 


85 


4 


2696 


42.0 


Sc 


Sy2 


85.8 


(A) 


24.3 


11 


11.38 


0.19 


N 3877 


Ilh46m07.80s 


+47d29m41.2s 


85 


5.1 


894 


17.7 


Sc 


HII 


40.7 


(A) 


5.19 


2.2 


10.74 


1.53 


N 4013 


Ilh58m31.38s 


+43d56m47.7s 


84 


4.7 


835 


16.0 


Sb 


HILLINER 


38.2 


(A) 


4.25 


1.8 


10.79 


1.34 


N 4096 


12h06m01.13s 


+47d28m42.4s 


82 


6.4 


577 


10.32 


SABc 


HII 


51.4 


(A) 


1.77 


0.75 


10.59 


0.40 


N 4157 


12hllm04.37s 


+50d29m04.8s 


90 


7 


771 


15.6 


SABb 


HII 


185 


(A) 


9.18 


3.8 


11.00 


1.19 


N 4192 


12hl3m48.29s 


+ 14d54m01.2s 


83 


8.7 


-142 


13.55 (V) 


SABb 


HII,Sy,LINER 


23 




3.3 


1.49 


11.15 


1.44 


N 4217 


12hl5m50.90s 


+47d05m30.4s 


78 


5.1 


1028 


20.6 


Sb 


HII 


120 


(A) 


11.8 


5.3 


11.03 


0.95 


N 4244 


12hl7m29.66s 


+37d48m25.6s 


90 


15.8 


247 


3.30 


Sc 


HII 


20.3 


(C) 


<0.03 


< 0.01 


10.02 


0.39 


N 4302 


12h21m42.48s 


+ 14d35m53.9s 


89 


4.7 


1118 


19.41 (V) 


Sc 


Sy.LINER 


34.8 


(A) 


5.08 


2.3 


10.90 


3.60 


N 4388 


12h25m46.75s 


+ 12d39m43.5s 


79 


5.6 


2607 


16.60 (V) 


Sb 


Sy2 


119 


(A) 


4.71 


2.6 


10.94 


1.56 


N 4438 


12h27m45.59s 


+ 13d00m31.8s 


78 


9.1 


259 


10.39 (V) 


Sa. 


LINER 


63.3 


(A) 


0.90 


0.32 


10.91 


2.67 


N 4565 


12h36m20.78s 


+25d59ml5.6s 


90 


16.2 


1228 


27.1 


Sbc 


Sy3,Syl.9 


134 


(D) 


15.7 


7.9 


11.80 


1.00 


N 4594 


12h39m59.43s 


-Ild37m23.0s 


79 


8.4 


1127 


12.7 


Sa 


LINER,Syl.9 


93.4 


(A) 


1.59 


0.68 


11.08 


0.32 


N 4631 


12h42m08.01s 


+32d32m29.4s 


85 


14.7 


613 


7.55 


SBcd 




772 


(C) 


8.49 


2.9 


10.74 


0.41 


N 4666 


12h45m08.59s 


-00d27m42.8s 


76 


4.2 


1516 


27.5 


SABc 


HII, LINER 


434 


(A) 


53.8 


20 


11.13 


0.54 


N 4845 


12h58m01.19s 


+01d34m33.0s 


81 


4.8 


1228 


16.98 (V) 


Sab 


HILLINER 


43.3 


(A) 


5.17 


1.8 


11.32 


0.49 


N 5084 


13h20ml6.92s 


-21d49m39.3s 


90 


12.5 


1728 


23.4 


SO 


poss. LINER 


45.9 


(A) 


0.49 


< 0.21 


11.91 


0.29 


N 5297 


13h46m23.68s 


+43d52m20.5s 


89 


5.3 


2404 


40.4 


SABc 




23.0 


(A) 


8.09 


3.9 


11.28 


0.79 


N 5775 


14h53m58.00s 


+03d32m40.1s 


84 


3.9 


1582 


28.9 


Sbc 




280 


(A) 


38.1 


14 


11.09 


0.67 


N 5792 


14h58m22.71s 


-01d05m27.9s 


81 


7.2 


1930 


31.7 


SBb 


HII 


51.8 


(A) 


17.2 


6.4 


11.47 


0.52 


N 5907 


15hl5m53.77s 


+56dl9m43.6s 


90 


11.2 


666 


15.26 


Sc 


HII 


50.5 


(A) 


5.56 


2.7 


11.32 


0.26 


U10288 


16hl4m24.80s 


-00dl2m27.1s 


90 


4.9 


2046 


34.1 


Sc 




26.1 


(A) 


2.55 


1.3 


11.07 


0.23 



Note. — Data are from[Tullv ( 1988j) unless otherwise indicated. 
a From the NASA Extragalactic Database (NED). 

inclination, i — 3deg + cos -1 (y / '((b / a) 2 — 0.2 2 )/(l — 0.2 2 )), where b/a is the minor to major axis ratio. 



c Observed blue diameter at the 25th mag/arcsec 2 isophotc. 
d Hcliocentric radial velocity. 

c Distance from NED, assu ming Hq — 73 km s~ 1 Mpc — 1 and correcting for Virgo Cluster and Great Attractor perturbations, except for Virgo Cluster galaxies 
which have been taken from (Solancs ct al. 2002); the latter arc indicated by (V). 

f From H YPERLEDA, [http : //leda . univ-lyonl ~frj\ see also lPaturel et al. I J2003f) . 

s Flux density at 1.4 GH z from refs: (A) NRAO VLA Sky Survey (NVSS [Condon et al."|fl99ll , (B) IIrwin. English Sorathia I J1999I) , (C) IStrong et al. I |2004h , 
CD HCondon et al. I J2002T) . 

h FIR luminosity from L F jr = 4tt D 2 F fir , wher e F FIR - 1.26 X 10~ 14 (2.58 Sqq + Sioo) W m~ 2 for S in Jy fsee lSanders fe Mir abel 1996), representing the 
flux in the 40.2 to 122.5 um range. Wc take S from Sanders ct al. (2003) or, if not available, from the most recent entry in NED, usually the IRAS Faint Source 
Catalog (FSC). N 4302 had no NED entry and was taken from the FSC directly. 

'SFR(Mq yr" 1 ; ) = 4. 5 X 10~ 44 L T /H(ergs s -1 ) JKennicutt Hl99Sf) where Ltir (3 to 1100 /J-m) is computed from S25, Sqq and Sioo according to the prescription 
in lDale fclleloujjl2002 ]) and t he flux densities are from the same sources as in (lie previous uote. These values may be refined in (he future with improved Hi data 
(see. e.g. lMurphv et al. Il2008h . 

^ (Log of) the total mass in units of Mq, where — W 2 d/(SG). W is the HI line width corrected for random gas motions and galaxy inclination, and d is the 

galaxy's linear diameter calculated from d25 and the galaxy's distance, corrected for the effects of projection and obscuration (values adjusted to our distance). 

k Density of galaxies brighter than -16 mag in the vicinity of the gala xy. Th e local density was determined on a 3D-grid at 0.5 Mpc spacing. No correction has 
been made for the fact that our distances differ somewhat from those of Tully (1988). 



Table 2 

Estimated Observing Time, rms & spatial coverage for CHANG-ES 



Array: i/o a 


Time/pt. b 


Thcor. rms c 


Confusion d 


®HPBW e 


Lin. Res. 


8las s 






(/ijy beam -1 ) 


(fijy beam -1 ) 


(arcscc) 


Min/max (pc/kpc) 


(arcmin) 


B: 1.5 GHz 


2 hrs 


7.7 




4.3 


69/0.88 


2.0 


C: 1.5 GHz 


30 min 


16 


9 


14 


224/2.85 


16.2 


6 GHz 


3 hrs 


2.4 


0.2 


3.5 


56/0.71 


4.0 


D: 1.5 GHz 


20 min 


20 


89 


40 


736/9.4 


16.2 


6 GHz 


40 min 


4.8 


2.3 


12 


192/2.4 


4.0 



NOTE. — These estimates represent a best-case scenario in which there is no RFI or other errors. 
a EVLA array configuration and central frequency. 

b On-source time per pointing, not including overheads. At 6 GHz, the field of view is smaller and the larger 
galaxies require more than one pointing. 

c Thcorctical noise, ignoring confusion, for each polarization product. Robust weighting ( Briggs 1995) is assumed. 
At 6 GH z the bandwidth, B — 2 GHz and at 1.5 GHz, B = 500 MHz. See the EVLA Observational Status Summary 
(OSS) at http://evlaguides . nrao . edu/ index . php?t itle=Category : Status 

d Estimatcd confusion limits (from the OSS). These confusion limits should only apply to total intensity, I, but 
not to polarization products, Q, U or V. 

c Spatial resolution, assuming uniform weighting. Actual values will be larger when robust weighting is used. 
Given the wide frequency bands and possibility of a variety of uv weightings, the range of spatial resolution is 
actually greater than noted here. 

f Minimum and maximum linear resolution for minimum and maximum galaxy distances, respectively, correspond- 
ing to the spatial resolutions of the previous column. 

8 Largest angular size scale detectable by the array. 
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